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Abstract 

A simple quark-diquark model for nucleon and A structure is used to 
calculate leading twist light-cone fragmentation functions for a quark to 
inclusively decay into P or A. The parameters of the model are deter- 
mined by fitting to the known deep-inelastic structure functions of the 
nucleon. When evolved from the initial to the final Q 2 scale, the calcu- 
lated fragmentation functions are in remarkable agreement (for z > 0.4) 
with those extracted from partially inclusive ep and e + e~ experiments at 
high energies. Predictions are made, using no additional parameters, for 
longitudinally and transversely polarized quarks to fragment into p or A. 
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Fragmentation of partons produced in a high-energy process into definite 
hadronic states is a problem in QCD of considerable current interest. Parton 
fragmentation, in a sense, is dual to the process of uncovering hadron structure 
using high momentum probes. This duality has recently been emphasized by Jaffe 
and Ji [1-2]. These authors have made a detailed exploration of the spin, chirality, 
and twist structures of fragmentation using the fact that parton fragmentation 
functions can be expressed in QCD as matrix elements of quark and gluon field 
operators at light-cone separations. Indeed, a one-to-one correspondence between 
fragmentation functions and parton distribution functions can be established. 
While this is a useful step towards understanding the nature of the hadronization 
process, it certainly does not solve the problem — the calculation of distribution 
functions belongs to the realm of non-perturbative QCD and therefore can only 
be modelled. Fragmentation functions are even harder to model theoretically, 
and existing models using strings and shower algorithms etc. [3] are complicated 
and involve many parameters. 

In this paper we explore an alternate, and much simpler, description for 
the conversion of a fast quark into specific hadrons. The starting point is an 
expression bilinear in quark field operators which defines the unpolarized and 
polarized fragmentation functions in terms of the light-cone momentum fraction 
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z = P + /k + . (seefig-1), 

4 M£l p + = f ^ e - iA/2 Tr < | 7+^(0) | PX >< PX | i>(\n) | > (1) 

4^Ap+S. n = f — e - lX ' z Tr < | 7 + 7 5 ^(0) | PS\\X >< PS\\X | $(An) | > 

(2) 

and, 

h ( z) f d\ 

A^^S ± p + = / — e-^ z Tr < | 7 + 7 ± 7 5 ^>(0) | PS ± X >< PS ± X | ^(An) | > 

Z J 2,71 

(3) 

In the above, and n M are null vectors with p 2 = n 2 = p~ = n + = and p.n = 1. 
In terms of these, the four momentum of the produced hadron, whose rest frame 
we shall take as our reference frame, is P M = p^ + |M 2 n M . The spin vector of the 
produced spin-± hadron is = S.np>* + S.pn^ + MS^, and 7 + = l/y/2 (7 + 7 s ). 
The gauge n.A = A + = is used. A summation over X is implicit and covers 
all possible states which can be populated by quark fragmentation. The quark 
operators are at equal x±_ and x + but are separated by a variable light-cone 
distance x~ = An~. To interpret fi(z) physically, define a composite operator 
C\P) which creates a hadron of a specific type and momentum from the vacuum, 
I P >= C\P) I >. Using completeness of X, it is easy to show that, 



*<«> ->* s- r 1 < ( :z m> « 
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where 



[dk] 



d 2 k T dk + 
(2tt) 3 2A;+ 



(5) 



This shows that fi(z) is the probability of finding a given hadron with fixed z in 



a quark, irrespective of the transverse quark momentum. 

The matrix elements in eqs.1,2 need to be modelled as they cannot be 
calculated ab-initio. Perhaps the simplest assumption is that the quark fragments 
into a baryon and anti-di quark (fig. lb) . This may be reasonable provided that z is 
not too far from 1, i.e. the produced hadron carries away most of the momentum 
of the quark. The amplitude for the process is, 



where <£> is the vertex that connects the quark with the two outgoing particles. 
$ is a matrix in Dirac space, and its most general form is rather complicated 
since it involves several unknown form factors. Instead, we shall be guided by 
the investigations of Meyer and Mulders[4], and Melnitchouk, Schreiber, and 
Thomas[5], who have calculated nucleon structure functions in the diquark model 
and obtained rather good fits to deep inelastic scattering at all data except small 
x. Following ref[4], we take 



< ps | i> | o >= u(ps)$ 



(6) 



^ — m 



$ = 0(P 2 ,p' 2 ,A; 2 ) J, 



(7) 
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where X is the unit matrix and only scalar diquarks are included in the vertex. 
Inserting eqs.6,7 into eqs. 1,2,3 and neglecting the quark mass m in the numerator 
yields, 

1 r oo dk 2 M 2 + z 2 k 2 , , |2 



and, 



fl{z) - MT^TjJo (2*) 2 {k 2 - m 2 ) 2 ^ ( 8 ) 



A . . 1 r°° dk\ M 2 — z 2 k\. , l2 



f , v M 2 r°° dk\ 1 . , |2 

hl( * ) = 8(1^)70 (10) 



The fragmenting quark, which is highly off-shell, has its k 2 given by, 



z 1 — z 



The integrals are logarithmically divergent and need to be damped by a suitable 
form factor. Following ref. [5] we choose to be, 



m = N W^y- (12) 



Here A is a mass parameter whose choice will be decided upon later. N is a 
normalization constant which is fixed by the requirement that the quark state be 
normalized covariantly, < k \ k' >= (2ir) 3 2k°5 3 (k — k'). 



Let us now concentrate on the spin-isospin structure of the produced baryon. 
In the notation of Meyer and Mulders [4], the SU(6) wavefunction for a spin- up 
proton is: 

Ip^^ifiS + ^gi^-^i^-ljT?-^! hi (13) 

S and T are quark combinations representing scalar and vector diquark respec- 

T T 

tively, e.g. T\ —uu, etc. There is, of course, no explicit coupling to vector 
diquarks in the vertex eq.7 although, at the expense of some complication, it 
could be put in. We follow instead the easy route suggested by the calculations 
of ref.[4] wherein it was assumed that the S(x) and T(x) distribution functions 
are of identical form but differ because the S and T diquarks have somewhat 
different masses m s and Tfly . CIS well as mass parameters A s and A„. The N — A 
mass difference leads to m v — m s ps 200 MeV. Following a similar logic, define 
S(z) and T(z) to be fragmentation functions where the emitted anti-diquark is 
repectively S and T, 

S(z) = f\{z) with m d = m s , A = A s 

T(z) = f\{z) with m d = m v , A = A v . (14) 

The polarized quantities AS(z), AT(z) are defined similarly with f\{z) replaced 

by <ji{z) and AiS(z), A]T(z) are defined with fi{z) replaced by h\{z) in the 
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above. The u and d fragmentation functions are readily seen to be, 



/£(*) = \s(z) + \f(z) 

/£(*) = \m 

h P M = l^S^-^fiz) 

h p ld {z) = -\^f{z) (15) 



The same approch can be followed for the production of a A hyperon from 
u, d, or s quarks. The hyperon wavef unction is, 

| A T>= -L= (u T°- d T° - V2 i T 1 + y/2 d T x + I S°+ dS° - 2 I S°^j . 

(16) 

Here S and T are diquarks with one s and either a u or d quark, while S is a ud 
system as before. Define new fragmentation functions, 

S(z) = fi{z) with m d = m s , A = A s 
T(z) = fi(z) with ma = m?, A = A? 
S(z) = fi(z) with m d = m s , A = A s . (17) 

The quantities AS(z), AT(z), AS(z) are defined similarly with fi(z) replaced 



by gi(z) and AiS(z), AiT(z), AiS(z) are defined with fi(z) replaced by hi(z). 
In terms of these, the fragmentation of u, d, and s into A is given by, 

/£(*)=/&(*) = fa) + 

ft(*)=ls(z), (18) 
and for the polarized quantities, 

{±Az) ~ Ait(z)) 
hUz)=l^iS(z) (19) 

We now turn to a discussion of numerical results and comparisons with 

fragmentation data, where available. The diquark model for fragmentation of u, d 

quarks to protons needs, as input, the scalar and vector diquark masses m s , m v as 

well as the cutoffs A s , A„. Consistent with the conclusions of Melnitchouk et.al[5], 

we find that a reasonable fit to nucleon structure functions can be achieved with 

the choice m s =900 MeV, m^=1100 MeV and A s =840 MeV, A„=925 MeV. This 

fixes all the parameters needed for the model, and one may use it for calculating 

fragmentation rates into the protons. Fig.2 shows ff u (z) and ff d {z) calculated 

using eqs.8,12. Since these are scale dependent quantities, one must specify the 
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scale as well. It appears reasonable to take the initial scale to be rrid + m p , where 
trip is the proton mass. This, of course, is the minimum off-shell mass of the 
fragmenting quark. Evolution of the fragmentation functions to the experimental 
scale can be performed exactly as for quark distribution functions. For example, 

/^//W) = Ejf J (^) (20) 

where P^j (x, /i) is the Altarelli-Parisi function for the splitting of the parton 
of type % into a parton of type j with longitudinal momentum fraction x. In 
principle, the sum in eq.20 extends over gluons and anti-quarks as well. At large 
z this is hopefully small and so we ignore this, including only P u -> u (x, //) and 
Pdr+d (x, fi) where, 

p„ m .2^4(^) + . (21) 

The + function is defined in usual way, 

f(x) + = f(x) - 5(1 - x) f 1 dx' f(x'). (22) 

Jo 

The evolved fragmentation function f£(z, Q 2 ) is compared in fig.3 with the EMC 
data [6] extracted from ji — p and /x — D scattering. The agreement is fairly good 
even down to rather small values of z, where model has no reason to be valid. 
Without changing parameters, predictions for gi and hi are shown in figs. 4, 5 
respectively. 
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The calculation of u, d, s quark fragmentation into A hyperons proceeds 
identically with the sole change of increasing the diquark mass by M\ — Mp = 
186MeV if the outgoing diquark contains a strange quark. The SU(6) wavefunc- 
tions then determine the relative magnitudes of fi, gi, and hi, which are plotted 
respectively in figs. 6-8. That the produced A carries the spin of the fragment- 
ing s-quark is apparent from fig. 7. The transverse fragmentation function hi for 
u, d, s — > A, if measured, would be a good tool for uncovering the transverse 
distribution function hi(x) [1,10] of the proton. While no experimental spin data 
exists, data from e + e~ collisions at large cm energies can be compared with pre- 
dictions of the present model provided we assume that only u, d quarks lead to 
proton production, and only u, d, s quarks lead to A production. Fig. 9 contains 
a comparison of theory to experiment. Again the evolution has been performed 
to a final scale of Q 2 = 80 GeV 2 . 

To conclude, we have investigated a simple fragmentation model for a quark 
to go into a P, A and the appropriate anti-diquark, using a vertex with a suitable 
form-factor determined by fitting to deep inelastic data. The good agreement 
with experiment— even down to rather small z values where this agreement is 
surely fortuitous— suggests that one is perhaps using the right effective degrees 
of freedom for P, A production from quarks. At the same time, it is obvious 
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that one is far from a complete description of fragmentation; the model gives 
zero chance for anti-proton production from a quark. Neverthelesss, it has good 
predictive power, and we have calculated the spin-fragmentation functions g\ and 
hi for P, A production. It will be interesting to see how well these are eventually 
borne out by experiment. 
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Figure Captions 

1. The fragmentation of a quark with momentum k (fig. la) into a specific 
hadron (P, A, • • •) is modelled (fig. lb) with a quark-diquark-hadron vertex. 

2. Fragmentation functions for u and d quarks to go into a proton, calculated 

in the diquark model, evolved from the initial scale Ql = (m^ + m p ) 2 to 
Q2 = 80 GeV 2 

3. Evolved diquark model fragmentation function f£(z, Q 2 ) compared against 
EMC [6] data extracted from \i - P and /i - D scattering. 

4. As in fig. 2, except that the fragmentation function §i is plotted here. 

5. As in fig. 2, except that the fragmentation function hi is plotted here. 

6. Diquark model calculation for u, d, s quarks to fragment into A hyperons. 
The initial scale Ql = (m d + m^) 2 and final scale Q 2 = 80 GeV 2 fragmen- 
tation functions are shown. 

7. As in fig. 6, except that the fragmentation function gi is plotted here. 

8. As in fig. 6, except that the fragmentation function hi is plotted here. 

9. Production of protons and hyperons in e + e~ collisions at cm energy of about 
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30 GeV. The solid curves are predictions of the diquark fragmentation 
model, and the data is from TPC[7], HRS[8] and MARK II[9] collaborations. 

Figures may be obtained by writing to 
hafsa%png-qau%sdnpk@sdnhq. undp . org 
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